Unscreened universality class for superconductors with columnar disorder 
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The phase transition in a model for vortex hnes in high temperature superconductors with colum- 
nar defects, i.e., linearly correlated quenched random disorder, is studied with finite size scaling and 
Monte Carlo simulations. Previous studies of critical properties have mainly focused on the limit of 
strongly screened vortex line interactions. Here the opposite limit of weak screening is considered. 
The simulation results provide evidence for a new universality class, with new critical exponents that 
differ from the case of strong screening of the vortex interaction. In particular, scaling is anisotropic 
and characterized by a nontrivial value of the anisotropy exponent = The exponents we 

find, ( = 1.25 ± 0.1,i^_L = 1.0 ± 0.1, z — 1.95 ± 0.1, are similar to certain experimental results for 
YBCO. 

PACS numbers: 74.60.-w (Type-II Sup.), 05.70.Fh (Phase Trans.), 75.40.Mg (Num. Simulations) 



Columnar defects have proven very effective at pin- 
ning vortex lines in high temperature superconductors 
(SC), resulting in significant increases in the critical mag- 
netic fields and currents 0. The columnar defects are 
often produced as damage tracks by irradiating the sam- 
ple with high-energy heavy ions. Such correlated disor- 
der strongly affects the phase transition from the SC to 
normal state and there is an ongoing theoretical and ex- 
perimental effort to investigate the phase diagram and 

§)perties of the transition 2, 3, 4, 5, 6, Z, ,8, 9, 10, 11, 
El El El El, E3 • Below a critical temperature the 
vortices freeze into a glassy SC state, known as a Bose 
glass, which replaces the Abrikosov vortex lattice of the 
clean system. Upon heating the system undergoes a con- 
tinuous phase transition to a resistive vortex liquid state. 
The transition is characterized by critical exponents, that 
have been measured in a number of experiments. The 
critical exponents are universal, i.e., they are common 
for all systems in the same universality class, and depend 
only on general features like the range of the interaction 
and distribution of the disorder. In spite of the fact that 
the screening length A in high temperature SC is quite 
large, most theoretical work on columnar defects has con- 
centrated on strongly screened systems with short range 
interactions [5 13. Il8| . In general, long range interactions 
often change the universality class of a transition, com- 
pared to short range interactions. It is therefore natural 
to study the universality class with long range interac- 
tions. In this paper we present a systematic study of the 
SC transition in systems with columnar disorder, an ap- 
plied magnetic field, and weak screening, and compare 
the results with experiments. 

The superconducting phase transition with colum- 
nar disorder is also interesting from the perspective 
of quantuin phase transitions in systems with disor- 
der H H ^1 0, I2I1. In fact, statistical mechanics 
of vortex lines in three-dimensions is closely analogous 
to world-lines of quantum bosons in (2-f 1) dimensions, 
where the vortex phase transition corresponds to the 



zero-tempcrature boson localization by substrate disor- 
der. The quantum bosons have static disorder in imag- 
inary time, corresponding to columnar disorder for the 
vortex lines. The quantum dynamical exponent ^ (relat- 
ing the diverging time and length scales at the transition) 
translates into an anisotropic scaling behavior for the vor- 
tex line problem, where correlations along the columnar 
defects diverge with a different rate than the perpendic- 
ular ones when approaching the Bose glass transition. 

The vortex phase transition with columnar disorder 
has previously been studied for the case of strongly 
screened vortex interactions There the correlation 

length exponent z/, anisotropy exponent and dynamic 
exponent z, take the values w 1, ^ = 2, z « 4.6 4J. The 
result C — 2 follows if the compressibility is nonsingular 
through the transition, and is hence believed to be ex- 
act. The effect of long range interactions in systems with 
correlated disorder has been considered in the context of 
the superconductor- insulator quantu m p hase transitions 
of disordered ultrathin films |l8[ lid . |2(j| . For long range 
planar 1/r-interactions, that act only between world line 
segments in the same xy-p\ane, and thus is completely 
local in the z-direction, the system is incompressible, 
with 1/ « 1 and C = 1 Ell- Early simulations [1 0| 
of more realistic models (from the vortex point of view) 
show evidence for a non-trivial value of the anisotropy ex- 
ponent close to 1. Some experiments ^1 0| have been 
interpreted as belonging to an incompressible universal- 
ity class with C w 1. 

The main task in this paper is to examine the super- 
conducting phase transition with columnar disorder in 
the limit where the screening length is very long, such 
that any effects of screening of the vortex interaction can 
be effectively neglected. In other words we set A 00, 
and thus examine a no-screening fixed point for the su- 
perconducting phase transition. The screening length in 
the HTS is typically A ^ 1000 A so for most of the phase 
diagram the assumption of no screening is reasonable. 
We stress that one possibility is that the true asymp- 
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totic critical behavior is effectively described by a fixed 
point with a screened interaction. However, even if this is 
the case, we expect that the scaling properties upon ap- 
proaching the transition is controlled by the unscreened 
model in a large region of the parameter space, before 
an eventual crossover to the true critical regime. In com- 
parison, this occurs in clean systems without an applied 
field |22| • Our model also assumes no fluctuations in the 
amplitude of the SC order parameter (London approx- 
imation), which is valid in a large region of the phase 
diagram well below Hc2 @- The interaction is taken to 
be a fully isotropic, three-dimensional (3D) long range 
interaction between vortex lines, that should apply for 
fairly isotropic systems. The form of the interaction 
makes our model different from most of the previously 
studied Boson-like models, with a planar interaction be- 
tween the vortex lines. We consider only the case when 
both the magnetic field and all the columnar defects are 
aligned in the z-direction. The anisotropy due to the 
field and the columnar disorder allows for the possibil- 
ity of an anisotropic correlation volume, where the cor- 
relation length in different directions can diverge with 
different exponents. In addition we consider dynamical 
aspects of the problem, which do not have a counter- 
part in the quantum boson problem, and compute the 
dynamical exponent z. This assumes that MC dynamics 
for vortex lines can be equated with real dynamics, which 
should apply close to the transition where the dynamics 
is slow and over-damped |23l|. Our values for Q, z differ 
from those of the previously studied models, which sug- 
gests that the present model belongs to a new, distinct 
universality class. 

Our starting point is the Ginzburg-Landau theory 
for the superconducting complex scalar order parameter 
field, ■(/'(r) = li/'l expi0(r). 
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where $o is the magnetic fiux quantum, B = V x A is the 
magnetic flux density and A is the magnetic vector po- 
tential, and H is the applied magnetic fleld. By fairly 
standard manipulations B this model can be trans- 
formed into a model involving only the vortex degrees 
of freedom. In the London approximation, i.e., neglect- 
ing fluctuations in the amplitude of '(/'(r), a discretized 
form of the vortex Hamiltonian reads 



^ = ^ E ^(r-r')q(r) •q(r') + i J2 D{v^)q.{vf 



(2) 



where the sums over r run through all sites on a simple 
cubic lattice with fl = L x L x sites, and periodic 
boundary conditions in all three directions. The vortex 



line variables are specified by an integer vector field q(r) , 
such that the fj, = x,y, z component is the vorticity on 
the link from the site r to r + . The partition function 
is Z = Tr exp(— _ff/T), where T is the temperature, and 
Tr denotes the sum over all possible integers q^, subject 
to the constraint that the discrete divergence V • q = 
on all sites, i.e., the vortex lines have no free ends. The 
vortex-vortex interaction is given by 
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V ^M(2-2cosfc^) + A-2' 



(3) 



where we choose units such that K — I, and consider the 
limit of no screening, corresponding to A = oo. Random- 
ness is included in the vortex model (0) in the second 
term in the form of a random core energy, which corre- 
sponds roughly to a random a(r) in the GL model. In 
real systems, each columnar defect gives approximately 
the same pinning energy, while they are located at ran- 
dom positions. On the lattice it is more practical to in- 
stead insert one column on each link of the lattice in the 
z-direction, but with a random pinning energy. We use 
a uniform distribution of random pinning energies in the 
interval < D{r±) < 0.8, where each D{r±) is constant 
in the z-direction to model correlated disorder. We model 
the net applied magnetic field as a fixed number of L'^/A 
vortex lines penetrating the system in the z-direction, 
corresponding to quarter filling, i.e. one vortex on every 
fourth link in the z-direction. To test for universality of 
our critical exponents, we also studied the case of L'^ /2 
vortex lines penetrating the system in z direction and 
disorder strength < D{r±) < 1.0. 

The MC trial moves consist of attempts to insert vor- 
tex line loops of random orientation on randomly selected 
plaquettes of the lattice. One MC sweep consists of one 
attempt on average to insert a loop on every plaquette, 
which we define as one MC time step. At — 1. The at- 
tempts are accepted with probability 1/(1 -|- exp AE/T), 
where AE is the energy change for inserting the loop. 
The initial vortex configuration is taken to be a regu- 
lar lattice of straight lines. To approach equilibrium we 
discard about 2 • 10^ MC sweeps (6.5 • 10** for the resis- 
tivity) before any measurements are taken, followed by 
equally many sweeps for collecting data. To verify that 
the warm-up time is long enough, we tried to vary the 
warm-up time between 5000 and 10^ sweeps for a single 
temperature, close to Tc, which gave no significant differ- 
ences in the final results. The results were averaged over 
up to 1000-2000 samples of the disorder potential. For 
the simulation of static quantities we use an exchange 
MC algorithm to speed up convergence Note that 
the exchange method can not be used for dynamic quanti- 
ties, since it involves large nonlocal moves in phase space. 
In the following we denote thermal averages by (...) and 
disorder averages by [...]. 

In the simulation, the helicity modulus and the RMS 
current is obtained by the following procedure |25|. An 
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extra term Hq = is included in the Hamiltonian, 

where is the total projected area of vortex loops added 
during the simulation |23- The helicity modulus in the 
direction fi (fi = x, z) is then given by 
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K 



(4) 



and the RMS current density is given by = 
where we use two different replicas in 
our simulations, denoted a and /3, to avoid any bias in 
the expectation values. We also calculate the linear re- 
sistivity p, in which case Hq is not included in H, by 
evaluating the Kubo formula 



t=-to 



(5) 



where t is MC time and oo, and the voltage is ~ 
AQ^. AQp is the net change in the projected vortex loop 
area during a sweep. In practice the summation time to 
is chosen large enough that the resistivity is independent 
of to- 
Next we consider the anisotropic finite size scaling re- 
lations used to extract critical properties from the MC 
data in Q . At the transition temperature Tc the corre- 
lation length in the xy-planes, ^, in the z-direction, 
and the correlation time, r, are assumed to diverge as 
^ - |r - Tc]^", 6 ~ T ^ The anisotropic 

finite size scaling ansatz 4] for the helicity modulus in z 
direction is 



(6) 



where d = 3 is the spatial dimensionality, and fz is a 
scaling function (all scaling functions will from now on 
be suppressed). Similarly, for the current density in the 
x-direction, we have 



2-d-C 



(7) 



The linear resistivity is given by p — E/ J, where E is 
the electric field and J the current density, and scales as 
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(8) 



In the simulations we consider a whole range of system 
sizes Lz for each L, since the anisotropy exponent C, that 
enters the aspect ratio, Lz ~ L'" is a priori unknown. 
For non- integral values of L^, we simulate two systems 
with nearest integer Lz and interpolate the results using 
linear interpolation. For the resistivity, both linear and 
logarithmic interpolations were tested and agree within 
error bars. Once we locate the correct value for C,, the 
argument Lz/L^ in the scaling functions can be made 
constant by setting Lz oc L^. 

We determine the critical temperature Tc for the phase 
transition and the critical exponents from MC data for 
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FIG. 1; MC data for the helicity modulus in the z-direction 
vs. temperature for different system sizes L. Tc ~ 0.0765 is 
estimated as the temperature where all curves intersect. Inset; 
Finite size scaling collapse of the data for fillings / = 1/2 
and / = 1/4, and for two disorder strengths (see text) using 
V = 1.0, C = 1-25. 
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FIG. 2: Finite size scaling collapse according to Eqs. Q and 
O of data for the helicity modulus in z-direction [main part] 
and RMS current in perpendicular directions [inset] at T = 
Tc The best collapse is obtained for T^ « 0.0765, C « 1-25. 



the helicity modulus in Eq. Q, by fits to the finite size 
scaling form in Eq. (jSJ. Figure ^ shows our estimate of 
Tc using the scaling form of the helicity modulus in the 
z-direction. By selecting the value for ^ that gives the 
best common intersection point, we obtain C, = 1.25 ±0.1 
and Tc = 0.0765, where the error bar on C is estimated 
by the interval outside which scaling gets considerably 
worse. The RMS current density in x-direction, J^, can 
also be used to calculate Tc- The result agrees, within 
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the error bars, with the result from Yz. 

Right at Tc it is possible to obtain a data collapse when 
plotting various quantities as a function of Lz/L'». Fig- 
ureElshow such plots for {L"^ /Lz)Yz and (L^L) Jmis- Ad- 
justing Tc and C, until the best collapse is obtained gives 
identical estimates, within the error bars, as the previ- 
ous method. This clearly demonstrates that the scaling 
is anisotropic. 

The inset of Fig. ^ shows a determination of the cor- 
relation length exponent v using finite size scaling, for 
data from two different fillings, / = 1/2 and / — 1/4, 
and two different disorder strengths. We use fits to Eq. 
(jHJ to obtain a data collapse for different system sizes 
(L = 4 - 8 for / = 1/2, L = 8 - 12 for / = 1/4) over 
an entire temperature interval around Tc. Both Tc and v 
were free parameters in the fit resulting in v — l.OitO.l, 
Tc = 0.0765 for / = 1/4, and T^ = 0.0747 for / = 1/2. 
The data collapse to a common universal scaling function 
indicates that the results are indeed universal. 
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FIG. 3: MC data for resistivity in the x and z-directions for 
different system sizes L at T = Tc. Straight lines are power 
law fits to the data points, with exponents z — C, = 0.65 and 
z-\-C, — 2 = 1.28. Inset: Finite size scaling plot of the resistivity 
as a function of the summation time to- 

Next we study the critical MC dynamics of the model. 
We calculate the resistivity p in x and z-directions from 
the Kubo formula in Eq. lP at T = = 0.0765. 
This gives us both a useful consistency test of the value 
C « 1.25 found above, and a value for the dynamic ex- 
ponent z. To verify that the summation time in Eq. 
© is long enough, we plot p vs. to in the inset of Fig. 
13 We also note that the correlation time is at most 
~ 10^ sweeps, which is much less than the equilibration 
times used in the simulations. Figure 13 shows the resis- 
tivity in the x and z-directions as a function of system 
size. We calculate C and z by making a power law fit 
of the data points in the figure to Eq. This gives 



TABLE L Summary of a few selected experiments and sim- 
ulations of high temperature superconductors with columnar 
disorder. 

Experiment/Simulation v z 

(K, Ba)Bi03 [101 1.1 ±0.1 ^2 5.3 ± 0.3 

TlaBaaCaCuaOs [8j 1.1 ± 0.2 1.9 ± 0.2 4.9 ± 0.2 

Bi2Sr2Cai-xYxCu208[llJ 1.04 ±0.06 ^ 2 5.28 ±0.05 

YBaaCusOr [14] 7^ l.Q ^ 1.1 ^ 2.2 

YBaaCusOy fl5] 0.9 ± 0.2 1.2 ±0.2 2.3 ± 0.3 

Simulations (A ^ 0) [4, fa 1.0 2 4.6 ± 0.3 

This work (A ^ oo) 1.0 ± 0.1 1.25 ± 0.1 1.95 ±0.1 



C — 1.3 ± 0.1, z = 1.95 ± 0.1, where the error bars are es- 
timated by the bootstrap method ^ in good agreement 
with C K, 1.25 found above. 

Finally we will compare our findings with some other 
results. The results for critical exponents, C = 1-25 ± 
0.1, 1' = 1.0±0.1,z= 1.95 ±0.1, imply that the linear re- 
sistivity scales as p ^ |T— Tel'* with Sx — v{z + 'i—d—C,) « 
0.7, Sz — v{z-\'l — d+C,) K, 1.3. For the nonlinear current- 
voltage characteristic at T = T^, we have E ^ ,P with 
Px - (l±z)/(d + C-2) « 1.3,p, = (C±z)/(d-l) « 1.6, 
for d — i. Table ^ shows a comparison between simu- 
lation results and results from a selection of transport 
experiments where the critical exponents studied in this 
paper have been measured. In the table we observe that 
several experiments agree very well with the previously 
known exponents for screened vortex interactions. No- 
tably, however, the table also shows that the results from 
two of the experiments on YBCO agree quite well with 
our exponents for unscreened 3D interactions, which sug- 
gests that they may effectively belong to the new uni- 
versality class considered here. Naively one may expect 
that the present results describe the transition when the 
screening length is much longer than the vortex spacing. 
For e.g. YBCO, where A ~ 1400 A, this corresponds to 
B > 0.1 T. The precise location of the crossover between 
weak and strong screening is unclear, as well as the pre- 
cise role of anisotropy and other parameters, and would 
be interesting to investigate further both experimentally 
and theoretically. 

In summary, we obtained scaling properties and nu- 
merical results for critical exponents that apply for the 
superconducting transition in systems with columnar de- 
fects in the limit of a long screening length, and com- 
pared with experiments. The critical exponents for this 
universality class differ considerably from the strongly 
screened case 13, B, and also from the case of planar 
1/r-interactions [ig. In particular, the new anisotropy 
exponent differs from the value C = 1 assigned to the uni- 
versality class of incom pres sible dirty bosons with planar 
long range interactions |18lll9ll20l |. Further work is moti- 
vated in order to further clarify the origin of the scaling 
properties obtained in different experiments, and when 
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the different models apply. Experimental studies to look 
for a crossover between unscreened and screened scaling 
behavior would also be interesting. 
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Goran Gustafsson foundation. 
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